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Abstract. We perform evolutionary calculations for a binary system with initial parameters: M 3 . 



1M© and 



M na ,i= 1.4 M© and P orbii (RLOF)= 1.27 d to produce observed binary parameters for the PSR J1740-5340. 
Our calculations support model proposed by D'Amico et al. (2001) in which this binary may be progenitor of a 
millisecond pulsar + helium white dwarf system. We propose observational test to verify this hypothesis. If the 
optical companion lack of carbon lines in its spectrum shows but the oxygen and nitrogen lines are present then 
our model correctly describes the evolutionary stage of PSR J1740-5340. 
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1. Introduction 

The millisecond pulsar PSR J1740-5340 was discovered 
during a systematic search of the Galactic globular clus- 
ter using the Parkers radio telescope (D'Amico et al. 2001 
a, b). The PSR J1740-5340 in the globular cluster NGC 
6397 is a member of a binary system with a relatively 
wide orbit of period 1.35 days, a companion with mass 
M c > 0.19 M Q , and is eclipsed for about 40% of its or- 
bit at a frequency of v=lA GHz. Spin-down age of PSR 
J1740-5340 {T sd =P pul /2P pul ) is -350 Myr and its surface 
magnetic field is B = 3.2 l0 19 (P pu iP P ui) 1/2 ~ 8 10 8 G. 
Ferraro et al. (2001) report the optical identification of 
the companion to the millisecond pulsar. They found that 
observed modulation is just — 0.2 mag; that can be repro- 
duced only if the companion has almost filled its Roche 
lobe and the orbital plane is nearly edge-on (i ~ 90°). 
These two requirements are accomplished by a compan- 
ion of mass in the range 0.19-0.22 M©, whose Roche lobe 
radius (1.32-1.42 R@) just matches the lower limit for ob- 
served radius R c ~ 1.3-1.8 i?©. D'Amico et al.(2001 a, b) 
and Ferraro et al. (2001) propose two hypothesis about 
the origin of this system. 

For explaining the large mass loss rate they proposed 
that the companion is an evolved star that spun up the 
pulsar to millisecond period. The mass accretion and spin- 
ning up of the pulsar would now be inhibited by the pulsar 
wind flux which could expel the matter overflowing from 



the Roche lobe of the companion. This would be the first 
confirmed example of a recently born millisecond pulsar. 

Another hypothesis suggested that the optical com- 
panion is a Main-Sequence (M-S) star perturbed by the 
energetic flux emitted from the millisecond pulsar. The 
only 5% of the impinging power would be needed to sus- 
tain the inferred mass loss rate M < 2 10 -11 M© yr~ x 
(D'Amico et al. 2001 b). This model should be more eas- 
ily applicable to a M-S star, as large convective envelope 
favours bloading of the star. 

Ergma, Sarna & Antipova (1998) investigated evo- 
lution and formation of short period binary millisecond 
pulsars near the so called bifurcation period (Pylyser & 
Savonije 1988). They found that if Roche lobe filling of 
the secondary occurs near bifurcation period then it is 
possible that the final orbital period (after detachment 
from Roche lobe) is near 1-2 days. 

In this research note we shall show that it is possible 
to produce observational parameters for the system PSR 
J 1740-5340 and that hypothesis proposed by D'Amico et 
al. (2001 a) really may work. If so we are observing binary 
right before formation of the detached He white dwarf + 
millisecond pulsar system. We also propose observational 
test to distinguish between the two hypothesis (M-S star 
vs. evolved star). 
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2. The evolutionary code 

The evolutionary sequences we have calculated are com- 
prised of two main phases: 

(i) detached evolution lasting until the companion fills its 
Roche lobe, 

(ii) semi-detached evolution. 



2.1. Detached phase 

The duration of the detached phase is somewhat uncer- 
tain; it may be determined either by the nuclear time-scale 
or by the much shorter time-scale of the orbital angular 
momentum loss owing to the magnetized stellar wind. 

2.2. Semi-detached phase 

In our calculations we assume that the semi-detached evo- 
lution of a binary system is non-conservative, i.e. the to- 
tal mass and angular momentum of the system are not 
conserved. The formalism which we have adopted is de- 
scribed in Muslimov & Sarna (1993, 1995). We introduce 
the parameter, /i, characterizing the loss of mass from the 
binary system and defined by the relations 



where L in is the intrinsic luminosity corresponding to 
the radiation flux coming from the stellar interior a is the 
Stefan- Boltzmann constant and R sg is the stellar radius. 
Lui is the millisecond pulsar radiation that heats the pho- 
tosphere, and is 



-'rot 



(4) 



and L rot is "rotational luminosity" of the neutron star 
due to magneto-dipole radiation (plus a wind of relativis- 
tic particles) 
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3c 3 ns 



2tt 



(5) 



where R ns is the neutron star radius, B is the value of 
the magnetic field strength at the neutron star and P pu i 
is the pulsar spin period. f 3 is a factor characterizing the 
efficiency of transformation of irradiation flux into thermal 
energy (in our case we take fs = 2 10~ 3 ). Note that in our 
calculations the effect of irradiation is formally treated by 
means of modification of the outer boundary condition, 
according to above the relation. 



M = M ss /i and M ns = -M sg (l - A), 



(1) 2.3. Pulsar spin evolution 



where M is the mass-loss rate from the system, M sg is 
the rate of mass-loss from the donor (secondary) star, and 
M ns is the accretion rate onto the neutron star (primary). 
The matter leaving the system will carry away its intrinsic 
angular momentum in agreement with formula 



J = M ns M 
J ~ h M so M 



yr 



(2) 



where M, M ns and M sg are the total mass of the sys- 
tem, the masses of the neutron star and donor star respec- 
tively. Here we have introduced the additional parameter, 
/ 2 , which describes the efficiency of the orbital angular 
momentum loss from the system due to a stellar wind 
(Tout & Hall 1991). In our calculations we assume /i=l 
(all transferred matter loss from the system) and /2=1 
(typical stellar wind - see Tout & Hall 1991). 

We also assume that the donor star, possessing a con- 
vective envelope, experiences magnetic braking (Mestel 
1968, Mestel & Spruit 1987), and as a consequence of this, 
the system loses its orbital angular momentum. 

In the semi-detached phase we have also included the 
effect of illumination of the donor star by the millisecond 
pulsar. In our calculations we assume that the illumina- 
tion of the component by the hard (X-ray and 7-ray) ra- 
diation from the pulsar leads to additional heating of its 
photosphere (Muslimov & Sarna 1993). The effective tem- 
perature, T e ff, of the companion during the illumination 
stage is determined from the relation 



Li, 



+ L M =4iraR 2 sg T* ff , 



(3) 



We also follow the rotational evolution of the neutron star 
which is determined mainly by the accretion rate and the 
temporal evolution of the magnetic field of the neutron 
star. For the temporal evolution of the magnetic field the 
model by Urpin & Muslimov (1992) has been used. 

As the neutron star spins up, the net accretion torque 
gradually decreases and the spin period reaches its equal- 
ibrium value given (see e.g. van den Heuvel 1977) by 



eg 
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When the spin period of the neutron star reaches P*, 
determined as 



P* = 1.9 B 1 ' 2 ( 



-1/4 



M, 



edd 



R 



1/8 



ms 



(7) 



where B 9 — B/10 9 G, M edd is the Eddington mass ac- 
cretion rate (for more details see Muslimov & Sarna, 
1993), then accreted matter is expected to be ejected from 
the system by the pressure of magneto-dipole radiation. 

2.4. Nuclear network 

Our nuclear reaction network is based on that of 
Kudryashov & Ergma (1980), who included the reactions 
of the CNO tri-cycle in their calculations of hydrogen and 
helium burning in the envelope of an accreting neutron 
star. We have included the reactions of the proton-proton 
(PP) chain. Hence we are able to follow the evolution of 
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the nuclei: 1 H, 3 He, 4 He, 7 Be, 12 C, 13 C, 13 N, 14 N, 15 N, 
uq 15q 5 16q, 17 and 17 F. We assume that the abun- 
dances of 18 and 20 Ne stay constant throughout the evo- 
lution. We use the reaction rates from: Fowler, Caughlan 
& Zimmerman (1967, 1975), Harris at al. (1983) and Pols 
et al. (1995). 

3. Results of calculations 

We perform calculations for the following system parame- 
ters: M sg4 =lM Q , M nSii =lA M Q , P orbii (RLOF)= 1.1-1.5 
d and chemical composition: X=0.7, Z=0.003. In Fig. 1 we 
show the mass-loss rate from the donor star versus sec- 
ondary mass. If Roche lobe filling of secondary occurs near 
bifurcation period (the bifurcation period separates the 
formation of the converging systems: Pf<P or b.i(RLOF) 
from the diverging systems: Pf>P or b i(RLOF), where Pf 
and P or b,i(RLOF) are the final and initial orbital periods, 
respectively) then as it was already found by Tutukov et 
al. (1985) the orbital evolution is quite different then in 
the case when the secondary Roche lobe filling occurs far 
from bifurcation period. At the beginning, mass loss rate 
is high and the pulsar will spin up to current spin period, 
then the mass loss rate quickly decreases. Near bifurcation 
period two time scales: nuclear time-scale of the secondary 
and the angular momentum loss timc-scalc arc very close 
each others. Therefore, during the binary evolution orbital 
period changes are insignificant in comparison with initial 
orbital period (Tutukov et al. 1985, Ergma et al. 1998). 

For various accretion rate we can estimate the opti- 
cal depth of the stellar wind due to free-free absorption 
(Illarionov & Sunyaev 1975) 
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Fig. 1. The evolution of the mass transfer rate versus sec- 
ondary mass. The horizontal line below which Tff is less 
than 1 is shown. To the right from vertical line the pulsar 
is in "propeller phase" 
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(log[Fe/H]) -2 and -3. The theoretical value of M„=4.44 
mag agree perfectly with Taylor at al. (2001) determina- 



In Fig 1. for secular mass loss rate the line below Tff is 
less than unity and the region where pulsar is in "propeller 
phase" (right from vertical line) are shown. Although the 
optical depth for secular mass loss is less than unity due to 
"propeller phase" plasma inside the system may be very 
clumsy with much higher optical depth and free-free ab- 
sorption alone could be viable eclipse mechanism. 

The model which provides the best fitting for the ob- 
served physical and orbital parameters is for P or b^=1.27 
d. During semi-detached evolution, near P orb =1.35 days 
secondary mass has decreased to M sg ~ 0.2 MQ, its 
radius is ~ 1.38P© and effective temperature 5620 K. 
Both values fit well with observed values (D'Amico et 
al. 2001 a, Ferraro et al. 2001). From Taylor et al. 
(2001) we know that 771^=16.66 mag for optical com- 
ponent of PSR J1740+5340. From Harris (1996) we 
find that distance modulus to NGC 6397 is 12.25 mag. 
These numbers gives absolut maginitude M V =AA1 mag. 
From our theoretical temperature and radius determi- 
nation we calculated M^oi — 4.166 mag. The bolomet- 
ric correction from Houdashclt, Bell & Sweigert (2000) 
gives BC=0.272±0.007 for assuming metalicity between 



tion. 

Fig. 2 a) and b) depict the temporal behaviour of the 
neutron star spin period, P pu i, as well as P eq and P*, in 
the case when B=5 10 8 and 8 10 8 G, respectively, at the 
beginning of the accretion stage. Fig. 2 shows that, prior 
to reaching its equilibrium value before the spin period 
reaches its equilibrium value (P e(? ), the spin period de- 
creases to the value of P* and the condition P pu i < P* is 
satisfied. After this, the evolution of the neutron star spin 
period is determined by magneto-dipole braking. When 
the spin period of the neutron star reaches P* (eq. 6) then 
the accreted matter is expected to be ejected from the 
system by the pressure of magneto-dipole radiation. It 
happens when M sg ~ 0.92 M , M sg ~ 6 10~ 9 M© yr^ 1 
for B=5 10 8 G and when M sg ~ 0.95 M , M sg ~ 
6.4 10~ 9 M Q yr- 1 for B=8 10 8 G, respectively. It should 
be pointed out that, in both cases the mass of the neutron 
star will grow, with ejection of accreted matter stopping 
at the value < 1.5 M Q at the beginning of the ejection 
phase. 

Note also, that for binary system parameters as above, 
pulsar spin period is 3.5 and 5 ms for 5 10 8 and 8 10 8 G, 
respectively. This will suggest that the value of 8 10 8 G 
given by D'Amico et al. (2001 b) is overestimated. 
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Fig. 2. The evolution of the spin period (P pu i) of the pul- 
sar, equilibrium period (P eq ) and P* versus the mass of 
the secondary for two values of magnetic field strength 
5 10 8 G (a) and 8 10 8 G . The horizontal dotted line 
shows the spin period of PSR J1740-53. The vertical lines 
show the mass range of the secondary star. 
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Fig. 3. The evolution of the red giant surface abundances 
of H, He, C, N and O as a function of the subgiant mass 
M sg . The vertical lines show the mass range of the sec- 
ondary star. 



4. Observational test 

We propose that observations of abundances of C, N, O el- 
ements in the spectra of optical companion of PSR J1740- 
5340 may give additional information about the nature 
of the secondary being either i) perturbed low-mass M-S 
star or ii) evolved star (helium) which has lost its envelope 
during the accretion phase. In Fig. 3 we present evolution 
of the optical companion surface chemical composition as 
a function of the secondary mass (for case ii)). Near or- 
bital period P or b= 1-35 d for the evolved star, carbon and 
oxygen are depleted by -2.6 dex and -0.8 dex (in compar- 
ison to cosmic abundances), respectively. Only nitrogen 
lines will show cosmic abundances (see Fig. 3). In the case 
of the low-mass M-S star C, N, O abundances will be as 
for the halo and disc metal poor ([Fe/H] < -1.5) dwarfs 
(Carretta, Graton & Sneden 2000 and references therein). 
In the case of the evolved star (helium) N will be over- 
abundant (approxymately solar) when O abundance will 
be as for metal poor dwarfs. It will be lack of carbon lines 
in the evolved star spectrum. 



5. Conclusions 

We accept D'Amico et al. (2001 a) hypothesis that the 
eclipsing binary millisecond pulsar PSR J 1740-5340 may 
be regarded as a recently born binary millisecond pulsar 
system. Our calculations shows that it is possible to build 
the evolutionary scenario for the binary system with the 
orbital parameters, secondary mass and radius which fit 
well those for PSR J1740-5340. We also suggest that the 
observed abundances of C, N, O elements in the spectra of 
optical companion of millisecond pulsar can tell us about 
evolutionary stage of secondary. If carbon is depleted then 
we really have a unique situation involving pre-Hc wd + 
millisecond pulsar system. 
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